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Abstract 

We present simple classical dynamical models to illustrate the idea 
of introducing a stochasticity with non-locality into the time variable. 
For stochasticity in time, these models include noise in the time vari- 
able but not in the "space" variable, which is opposite to the normal 
description of stochastic dynamics. Similarly with respect to non- 
locality, we discuss delayed and predictive dynamics which involve 
two points separated on the time axis. With certain combinations of 
fluctuations and non-locality in time, we observe a "resonance" effect. 
This is an effect similar to stochastic resonance, which has been dis- 
cussed within the normal context of stochastic dynamics, but with 
different mechanisms. We discuss how these models may be devel- 
oped to fit a broader context of generalized dynamical systems where 
fluctuations and non-locality are present in both space and time. 
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1 Introduction 



"Time" is a concept that has drawn a lot of attention from thinkers in virtu- 
ally all disciplines [1] . In particular, our ordinary perception is that space and 
time are not the same, and this difference appears in various contemplations 
of nature. It appears to be the main reason for the theory of relativity, which 
has conceptually brought space and time closer to receiving equal treatment, 
and it continues to fascinate and attract thinkers from diverse fields. More- 
over, issues such as the "direction" or the "arrow" of time [2] and complex 
time [3] are current research interests. 

It seems that there are other manifestations of this difference. One is the 
treatment of noise or fluctuations in dynamical systems. Time in dynamical 
systems, whether they are classical, quantum, or relativistic, is commonly 
viewed as not having stochastic characteristics. In stochastic dynamical the- 
ories, we associate noise and fluctuations with only "space" variables, such 
as the position of a particle, and not with the time variable. In quantum 
mechanics, the concept of time fluctuation is embodied in the time-energy 
uncertainty principle. However, time is not treated as a dynamical quantum 
observable, and clear understanding of the time-energy uncertainty has yet 
to be found [I]. 

Another difference seems to show up in our cognition of non-locality in 
space and time. Non-local effects in space are incorporated in physical the- 
ories describing wave propagation, fields, and so on. In quantum mechanics, 
the issue of spatial non-locality is more intricate, constituting the backbone 
of such quantum effects as the Einstein-Podolsky- Rosen paradox [5]. With 
respect to time, there have been investigations of memory effects in dynami- 
cal equations. However, less attention has been paid to non-locality in time, 
and behaviors associated with non-locality in time, such as delay differential 
equations, are not yet fully understood [6l [7J [8j [9] . 

Against this background, the main topic of this paper is to consider 
stochasticity and non-locality of time in classical dynamics through a pre- 
sentation of simple models. We discuss delayed and predictive dynamics as 
examples of non-locality in time. For stochastitiy, we present a delayed dy- 
namical model with fluctuating time, or stochastic time. We shall see that 
this combination of stochasticity and non-locality in time can exhibit be- 
haviors which are similar to stochastic resonance [TQl CHI EE2] , which arises 
through a combination of oscillating behavior and "spatial" noise and has 
been studied in variety of fields [13], [HI [HI HH [17] . 
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2 Delayed and Predictive Dynamics 



We start with a consideration of non-locality of time in classical dynamical 
models. The general differential equation of the class of dynamics we discuss 
here is as follows. 

^ = /(*(*), s(t)). (1) 

Here, x is a dynamical variable of time t, and / is the "dynamical func- 
tion" governing the dynamics. Its difference from the normal dynamical 
equation appears in t, which can be either in the past or the future, and 
t 7^ t in general. In other words, the change in x(t) is governed by /, not 
its "current" state x(t), but its state x at t. We can define t and x, as well 
as the function /, in a variety of ways. In the following, we will present two 
cases: delayed and predictive dynamics. 

Delayed dynamics can be obtained from the general definition by 

t = t-r, x(t)=x(t-r). (2) 

Here, r > is the delay, and the dynamics depend on two points on the 
time axis separated by r. Delayed dynamical equations have been studied 
for various applications [6j [7J [H [9] . 

Predictive dynamics, on the other hand, have recently been proposed 
[T8l [19] and, they take t in the future, i.e., t — t + rj. We call r\ > an 
"advance" . We also need to define the state of the dynamical variable x at 
this future point in time. Here, we estimate x such that 

x(t = t + rj) = v ( ^- + x{t). (3) 

This prediction is termed "fixed rate prediction" . Namely, we estimate x as 
the value that would be obtained if the current rate of change extends for 
a duration from the present point to the future point. Qualitatively, this is 
one of the most commonly used methods for estimating population, national 
debt, and so on. 

We also note that there are studies of equations called "functional differ- 
ential equations of the advanced type", or "advanced functional differential 
equations" j20l [21], [22] . They also are differential equations with advanced 
arguments, and we can obtain equations of this class from our general defi- 
nition by setting, 

x(i=t + 7]) = x{t + r]), (4) 
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with suitably chosen boundary conditions. The predictive dynamical equa- 
tions differ from this class of equations, as we allow flexibility in defining x 
based on a prediction scheme. 

We shall investigate the properties of these delayed and predictive dy- 
namical models through computer simulations. To avoid ambiguity and for 
simplicity, we will study time-discretized map dynamical models, which incor- 
porate the above-mentioned general properties of the delayed and predictive 
dynamical equations. 



x n+l = (1 - a)x n + f[Xn] 



(5) 



Here, a is a parameter controlling the rate of change. For a delayed map 
with a delay r, we define 

J'n = ■En—T) (6) 

while for the predictive map with an advance rj, we have 

Xfi = x n -\- Tj{x n X n — i). (7) 

We choose the Mackey-Glass function as the dynamical function (Fig. 1), 

fix 



i.e. 



/(*) 



(8) 



1 + X s 

where (3 and s are parameters. This function was first proposed for modeling 
the cell reproduction process and is known to induce chaotic behavior with 
a large delay [5J. 




Figure 1: Mackey-Glass function f(x) with f3 = 0.8 and s = 10. The straight 
line has a slope of a = 0.1. 
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Figure 2 shows examples of computer simulations of the delayed and 
predictive cases. The parameters are set so that without a delay or advance, 
t = 7] = 0, the model monotonically approaches the stable fixed point. The 
stability of the fixed point is lost as r, or 77, increases, giving rise to complex 
dynamics. Thus, non-locality in time can induce a complex behavior in 
otherwise simple dynamical systems. 

Now, we would like to make a few remarks. First, in the case of delayed 
dynamics, we need to decide on the initial function and delay. Analogously, in 
predictive dynamics, the prediction scheme and advance need to be specified. 
Common to delayed and predictive dynamical systems, both factors affect the 
nature of the dynamics. 

In addition, we can use linear stability analysis on both the delayed and 
predictive cases. This analysis can give an estimate of the critical delay or 
advance at which the stability of the fixed point is lost. However, the nature 
of the dynamics beyond these critical points is not yet clearly understood. 

For the case of delayed dynamics, with the addition of a suitable "strength" 
of noise, a behavior similar to stochastic resonance has been obtained [23J. 
This phenomenon, called "delayed stochastic resonance," has a different 
mechanism in the sense that it does not require external oscillatory signals or 
forces, but instead it uses a delay as the source of the oscillation to be com- 
bined with noise. An analogous resonance phenomenon has been observed in 
predictive dynamics with added noise, which is termed "predictive stochastic 
resonance" [19] . 



3 Stochastic Time 

We now turn our attention to the fluctuation of time, which we term "stochas- 
tic time" in the context of classical dynamical systems. As in non-locality, 
there are various ways to bring in stochasticity. We have found that stochas- 
tic time combined with delayed dynamics leads to phenomena similar to 
stochastic resonance. 

The general differential equation of the class of delayed dynamics with 
stochastic time is given as 

^ = f(x(t),x(i-r)). (9) 

Here, as in the previous section, x is the dynamical variable of time t, 
and / is the "dynamical function" governing the dynamics, r is the delay. 
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Figure 2: Examples of delayed (left column) and predictive (right column) 
dynamics for the Mackey-Glass map with a = 0.5, (3 = 0.8, and n = 10. For 
delayed dynamics, the initial condition is fixed at xq = 0.5 for the interval 
of (— r, 0). For predictive dynamics, the initial condition is Xq = 0.5 and 
Xi — (1 — h)xq + f(x ). The values of delay and advance r = r\ are (a) 0, (b) 
2, (c) 8, (d) 10, and (e) 20. 



6 



The difference from the normal delayed dynamical equation appears in i, 
which now contains stochastic characteristics, and these can be introduced 
in various ways. We will again focus on the following dynamical map system 
incorporating the basic ideas of the general definition given above. 

x n k+1 — f(x nk ,Xn k - T ), 

n k+1 = n k + £ k (10) 

Here, is the stochastic variable which can take either +1 or —1 with certain 
probabilities. We associate "time" with an integral variable n. The dynamics 
progress by incrementing integer k, and n occasionally "goes back" a unit 
with the occurrence of £ = —1. Let the probability of = — 1 be p for 
all k, and we set n = 0. Then, with p = 0, this map naturally reduces to 
a normal delayed map with n k = k. We update the variable x n with the 
larger k. Hence, x n in the "past" could be "re-written" as n decreases with 
probability p. 

We can make an analogy of this model with a tele-typewriter or a tape- 
recorder, which occasionally moves back on a tape. Figure 3 gives a schematic 
view. Based on the values of x n and x n - T , the recording device writes on 
the tape the values of x at a step, and "time" is associated with positions on 
the tape. When there is no fluctuation (p — 0), the head moves only in one 
direction on the tape and it records values of x for a normal delayed dynamics. 
With probability < p, it moves back a unit of "time" to overwrite the value 
of x. The question is how the recorded patterns of x on the tape are affected 
as we change p. 
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Figure 3: Schematic view of the model. 



We will keep the Mackey-Glass function as the dynamical function, and 
the map model becomes 



x 



v l oi^x nk -\- f{x nk — T ), 



n k+ i = n k + i 



fcj 
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/(*) 



/3x 



(11) 



1 + X s ' 



where a, (3, and s are parameters. With both a < (3 positive, and no 
stochasticity in time, this map has a stable fixed point with no delay. Linear 
stability analysis around the fixed point gives the critical delay r c , at which 
the stability of the fixed point is lost. 



A larger delay gives an oscillatory dynamical path. We have found, 
through computer simulations, that an interesting behavior arises when the 
delay is smaller than this critical delay. The tuned noise in the time flow gives 
the system a tendency for oscillatory behavior. In other words, by adjusting 
the value of p controlling £, one induces oscillatory dynamical paths. Some 
examples are shown in Figure 4. (The critical delay is r c ~ 22.5 with the set 
of parameters.) With increasing probability for a time flow to reverse, i.e., 
with increasing p, we observe oscillatory behavior in the sample dynamical 
path as well as in the corresponding power spectrum. However, as p increases 
further, the oscillatory behavior begins to deteriorate. In order to see this, 
we compute the "signal-to-noise" (S/N) ratio by using the ratio of the peak 
height to the background in the spectrum. Figure 5 illustrates this change 
in S/N, which reaches a maximum at an appropriately "tuned" value of p. 

Again, we see a phenomenon which resembles stochastic resonance. A 
theoretical analysis of the mechanism of our model is yet to be done. In 
particular, we note that the model has an intricate mixture of time scales 
of delay, the oscillation period of x n , and stochastic time. These factors are 
likely to be involved in the time scale analysis [24], but we leave consideration 
of them for the future. On the other hand, this resonance with stochastic 
time is clearly of a different type and new. We discussed only the Mackey- 
Glass function for consistency, but the same behavior is also observed in other 
delayed dynamical systems, such as ones with a negative feedback function 



We have also studied predictive dynamics with stochastic time, but so far, 
have not found behaviors similar to this resonance with delayed dynamics. 
Yet, the example above with delayed dynamics indicates that a combination 
of stochasticity and non-locality in time may lead to entirely new phenomena. 




(12) 



[25] 
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Figure 4: Dynamics (left) and power spectrum (middle) of delayed dynamical 
model with stochastic time. (The right column plots the values of stochastic 
time function of k.) This is an example of dynamics and associated 

power spectrum simulated with the model of Eq. (10) with a variable prob- 
ability of stochastic time flow p. The parameters are a = 0.03, (3 = 0.05, 
s = 10, and r = 15, and the stochastic time flow parameter p is set to (a) 
p = 0, (b) p = 0.15, (c) p = 0.3, (d) p = 0.4, and (e) p = 0.45. We used 
the initial conditions x n = 0.5(n < 0), and n = 0. The simulation had 
k = 10240 steps, and the values of x n for < n < L with L = 1024 were 
recorded at that point. Fifty averages were taken for the power spectrum of 
this recorded x n . The unit of frequency A is set as -j-, and the power -P(A) is 
in arbitrary units. 



Figure 5: Signal-to-noise ratio S/N at the peak as a function of the prob- 
ability of stochastic time flow p. The parameter settings are the same as in 
Figure 4 with (a) r = 10, (b) r = 15, and (c) r = 20. 

4 Discussion 

We could extend our model so that we have a picture of dynamical systems 
with stochasticity and non-locality on the time and space axes. The analyt- 
ical framework and tools for such descriptions need to be developed, along 
with a search for appropriate applications. 

An example of an appropriate application of temporal non-locality is mod- 
eling a stick balancing on a human fingertip. Recent experiments have found 
that most of the observed corrective motions occur on shorter time scales 
than that of the human reaction time [2H1 EE]- This may be the result 
of intricate mixtures of physiological delays, predictions, and physical fluc- 
tuations. Models incorporating special fluctuations and delays have been 
considered, but none have tried to include the effect of prediction. 

Another direction of development might be to extend the path integral 
formalism to allow stochastic time paths. The question of whether this ex- 
tension bridges to quantum mechanics and/or leads to an alternative under- 
standing of such properties as the time-energy uncertainty relations requires 
further investigation. 

Finally, if these models can capture some aspects of reality, particularly 
with respect to temporal stochasticity, this resonance may be used as an 
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experimental indication for probing fluctuations or stochasticity in time. We 
have previously proposed "delayed stochastic resonance" [23] , a resonance 
that occurs through the interplay of noise and a delay. It was theoretically 
extended [29], and recently, it was experimentally observed in a solid-sate 
laser system with a feedback loop |30j. We leave it for the future to see if an 
analogous experimental test could be developed with respect to stochasticity 
and non-locality of time. 
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